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1. Introduction 


Since the introduction of the Erdos-Renyi model [MlHS], random graphs have been studied extensively as 
central objects in probability theory [nun], computer science [UMIIST], and biology [T]I5]I32I]33) . In addition 
to furthering our understanding in these fields, random graphs arise naturally as models that reproduce many 
properties found in real-world networks [5]ll0lll2l[T^27ll38ll39j . In many random graph models, the degrees of 
the vertices concentrate around a single value with small fluctuations. However, in some real-world networks 
the behavior of the degrees is evidently different [2]. A graph is called scale-free if the number of vertices 
with degree j is proportional to for some positive /3. 

The spectral properties of random graphs are of great interest in applications as they relate to many 
combinatorial properties of the graph, to the mixing times of Markov chains, etc 0111]. In this paper, when 
we refer to any spectral property of a graph, we mean the eigenvalues and eigenvectors of its adjacency matrix. 
Many different statistics are used to study the spectra of random graphs. The statistics can be divided into 
two main types: bulk and edge. Bulk statistics involve the eigenvalues (and their respective eigenvectors) in 
the interior of the spectrum. More precisely, if we consider the empirical spectral distribution (esd) of an 
N X N random matrix, where a point mass of weight I/TV is placed at the location of each eigenvalue, then 
the bulk of the spectrum contains any interval where the limit of the esd has a density bounded away from 
zero. Edge statistics concern the behavior of the extremal eigenvalues at the edges of the spectrum. At an 
edge, the limiting density of the esd falls off to zero. 

When the esd converges in distribution, we call the limit a global law as this tells us that the density 
of eigenvalues in a macroscopic interval of length order one for a large hnite matrix is well approximated 
by the global law. At microscopic scales, we can ask about the density of eigenvalues in intervals of length 
order for positive e. A local law proves that the density of eigenvalue in these microscopic windows 

is also well approximated by the limiting distribution. We cannot expect this to hold for intervals of length 
less than order N~^, as the number of eigenvalues in such an interval does not grow with N, and so 
is optimal.To study this convergence of the density of the eigenvalues, the Stieltjes transformation, which 
is defined as := f (t — for a measure /x, is often used as this is equivalent to convergence 

in the esd. The parameter z = E + irj e C+ tracks the density of the eigenvalues at energy E for spectral 
windows of length rj. Sometimes the limiting distribution has a closed form, like the semicircle distribution, 
but for more complicated models, it is often specified as a functional self-consistent equation of the Stieltjes 
transformation. The bounds on the rate of convergence of the Stieltjes transformation further divide local 
laws. A strong law is optimal and provides the bound whereas weaker laws bound the convergence 

with a smaller power oi Nrj. 

While the global and local laws are model dependent, the fluctuations of the eigenvalues about the loca¬ 
tions predicted by the global law show substantial universality. In particular, the behavior of the gaps between 
eigenvalues and the n-point correlation function of the eigenvalues from many real-symmetric random matrix 
ensembles are often the same as the Gaussian orthogonal ensemble. This observation is formalized by the 
Wigner-Dyson-Gaudin-Mehta conjecture, or bulk universality conjecture, that states that the local statistics 
of Wigner matrices are universal in the sense that they depend only on the symmetry class of the matrix, 
but are otherwise independent of the details of the distribution of the matrix entries. This conjecture for all 
symmetry classes has been established in a series of papers [niiisHis] After this work began, parallel results 
were obtained in certain cases in [35ll36] . 

Edge statistics often exhibit universality too. The extremal eigenvalues are studied as a point process 
[7] and typical results involve the almost sure location that the largest eigenvalue converges to and the 
fluctuations about this location. Often these fluctuations, which are independent of the details of the model 
and depend only on the symmetry class, fall into the Tracy-Widom universality class. 
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Returning to random graphs, for the original Erdos-Renyi model, which is distinguished from Wigner 
matrices by its sparsity and discreetness, much has been proved. The degree distribution for the Erdos- 
Renyi model concentrates about Np with small fluctuations (so long as p is not too small) and the global 
distribution is given by the Wigner semicircle. A strong local law at the optimal scale for the density of the 
eigenvalues is known, so long as pN » (log A^), along with universality of both the bulk and edge statistics, 
when pN » [T7l[T8]. Recently, universality was proved to hold down to the scale pN » N'^ in . 

However, when the degree distribution is given by a power-law, the global distribution of the eigenvalues 
does not follow the Wigner semicircle—instead it follows a power-law [34]. This can be seen empirically in 
real-world data and in various models that produce scale-free random graphs. To state more precise results, 
we must focus on a model: in the Chung-Lu random graph model [EKTS], the expected degree sequence of 
the graph (di ,... ,(1^) on N vertices is specified and an undirected edge is added between vertices i and j 
independently with probability pij := didj/D^ where D := '^f.dk- Many different degree distributions can 
be produced by the Chung-Lu model, including power-laws with an exponent /3 > 3. For general degree 
sequences, the almost sure location of the largest eigenvalue is known [15] and the global law of the spectrum 
is obtained as the solution of a functional equation for its Stieltjes transform (see (12.131) 1 , which was derived 
in [33] using the replica method. 

The distinguishing feature of the Chung-Lu model is that the variance matrix has low-rank. Motivated 
by this (see Section Ej), we study the following more general model, where H = ( hij ) is a real symmetric 
N X N random matrix. The entries of H are independent (up to the symmetry constraint), have mean zero, 
and variance matrix S = (s^) with low-rank. The matrix H can also be sparse and we parameterize the 
sparseness with q = N'^ for k g (0,1], which we point out is different from the parameter in [T7l[T8]. The 
moments of the entries decay as ]E|/iy | ^ for fc > 3. Just like the Chung-Lu model, the 

global law of the spectrum is obtained as the solution of a functional equation for its Stieltjes transform 
(see Section (5]). We derive this self-consistent system in Section |31 In Section [S] we prove existence and 
uniqueness of the solution using the Brouwer fixed-point theorem and an elementary argument. Additionally, 
we prove the regularity and stability of the solution and provide its asymptotics. 

For this ensemble, we prove a local law in the bulk at the optimal scale rj » The result is parallel 

to that in m- The local law enables us to prove bulk universality, using a recent result in [3T]. Our bulk 
universality result is parallel to that in [50], and is proved using the same technique. We point out that in 
the non-sparse case, that is q = N, our results are contained in the results by Ajanki, Erdos, and Kruger [4]. 
The novelties in our paper are: we introduce sparsity into the ensemble; we analyze a different equation 
(12.131) than the quadratic vector equation (qve) introduced in [3], in order to handle possible singularities 
in Sy. 

In the appendix, we give a new proof to the existence and uniqueness of solution the QVE. We also study 
a general sparse ensemble, where we no longer require Sy to be low-ranked, but instead assume it to be flat 
in the sense that Sij ^ C/N. We state and sketch the proof of a local law and bulk universality of this 
ensemble. As mentioned above, Ajanki, Erdos and Kruger [4] have earlier proved the local law and bulk 
universality in the non-sparse setting {q = N). 

The layout of the paper is as follows: In SectionjS] we define the random-sign model, introduce some basic 
definitions, and precisely state our main results. In Section jS] we introduce some specific ways to produce 
degree sequences from a distribution tt that satisfy our assumptions. In SectionjT] we prove Theorem l2.6l the 
local law. Section |S] is devoted to proving the existence, uniqueness, and regularity of the solution to (12.131) . 
In Section |6l we use a Green’s function comparison argument and a result of [31] to prove the universality 
in the bulk for q » N’^. In Section jT] we consider a different model, where the entries only take the values 0 
and 1. We sketch the proof for the same local law and bulk universality for this model under the assumption 
that q « Appendix lAl is devoted to proving the existence, uniqueness, and regularity of the solution 
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to a general self-consistent system. In Appendix |B1 we state the local law and bulk universality for a general 
matrix ensemble whose limiting density is given by the self-consistent equations in Section lAl 


2. Definition and main results 


In the Chung-Lu random graph model [HIE], the expected degree sequence of the graph (di,..., d^) on N 
vertices is specified and an undirected edge is added between vertices i and j independently with probability 


P^J ■■ = 


didj 

~1T 


( 2 . 1 ) 


where D := Note that the probabilities ptj are only defined if max^ ^ D. Obviously, the exact 

degree of a vertex i is random, but its expected degree is exactly di. Under a particular choice of the expected 
degree sequence, which we specify in Section [31 the model can produce a scale-free graph for any fd > 3. 
Moreover, the original Erdos-Renyi model is recovered as a special case when the probabilities pij do not 
depend on i or j. In fact, under some mild assumptions on the expected degree sequence, a wide class of 
degree distributions can be produced by this model. Often the most interesting choices for pij lead to sparse 
graphs. 

To study the spectral properties of the adjacency matrices, denoted by A^v, of random graphs from the 
Chung-Lu model, we must center and rescale their entries; this produces a symmetric N x N matrix model 
Hn = The spectral properties cannot be studied if the graph is too sparse. We parameterize 

the sparseness with q = N'^ := D/N for k g (0,1], which we point out is different from the parameter 
in HZlllH] —in fact, gHERE = gTHERE- The parameter q is the order of the average degree of the graph and note 
that we require that k is positive. There are two ways to transform the entries to have expectation zero. 
The entries can be centered by subtracting their mean. Alternatively, each entry can be multiplied by an 
independent random sign (up to the symmetry constraint). Next, to ensure the density of the eigenvalues 
of converges to a distribution of order one, we rescaled by 

Centering and rescaling with the transformation hij = (uij —Pij)l-s/q leads to independent entries (up to 
the symmetry constraint) that have the following distribution 


P [hij = (1 - Pij) /sfq ] = p^J and P = -pijj^Jq ] = 1 - pij . (2.2) 

The random sign approach with the transformation hij := SijaijI^Jq (where Sij are independent random 
signs) again leads to independent entries, but they have the distribution 

P [hij = ±l/^/q ] = Py/2 and P [hij = 0] = 1 - pij . (2.3) 


The first matrix model has the variance matrix 


S' = 


didj {D - didj) 
qD^ 


which has rank two, and the second has the variance matrix 


S = 



which has rank one. 


(2.4) 


(2.5) 
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Thus, we are lead to the following generalized ensemble, which includes the above matrix models as 
examples. Denote Nat := {1,..., N}. We introduce an ensemble oi N x N symmetric random matrices 
The parameter N is large and we often omit explicitly indicating the dependence of 
various quantities on N. As discussed in Section [ 1 ] we motivate the ensemble by considering the adjacency 
matrices of random graphs from the Chung-Lu model, where edges are added independently between two of 
the N vertices, i and j, with varying probabilities pij (see Equation (l2.1D b 
Let Sij = E\hijf. We assume that 


Ehij = 0 and 


Nqk/2 


-1 


( 2 . 6 ) 


for fixed fc G N. We assume that the variance matrix S = is of low-rank and has the form 


1 V* (^) (^) 
Zj a 7) 


'j 


(2.7) 


fc=i 


where r is a fixed positive integer and ( 7 , 


(fc) 


)feeNr.ieNiv Satisfy the following assumptions: 


(i) 7= 1 for all k eNr and i G Nat; 

(ii) there is an M > 0 (not dependent on N) such that Zfe i(7i^^)^ ^ ^ ■ 


We also require that 


qstj < 1 , 


( 2 . 8 ) 


so that the matrix can be realized through the adjacency matrix of a random graph. 

Remark 2.1. It is possible that max^j Sij ^ 00 as ^ 00 , as shown in Subsection 13.21 In the non-sparse 
case q = N, this is impossible, since qsij ^ 1 implies ^ 1/A^. 

We denote 

(2.9) 

k 


Definition 2.2. If H = (/iij)i,jeNN satisfies the conditions above, we call H a generalized Chung-Lu ensem¬ 
ble. 


We are going to analyze the behavior of the resolvent of H, that is. 


G{z) = {H-z)-^ 


( 2 . 10 ) 


for z = A + i 77 G C+. 

We define the Stieltjes transform of a measure p as mp{z) '■= Denote the eigenvalues of II, 

ordered in increasing size, by (Ai, ... ,Xn). We are interested in the asymptotic spectral statistics of H, so 
we consider the empirical spectral distribution of H, defined as /iAi •= Yji and its Stieltjes transform 


mN{z) 


Aijv(dcc) _ J_ y> 1 
X — z N ^ \i — z 

I 


( 2 . 11 ) 
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In terms of the resolvent (|2.1Qp . tun can be written 


1 

«T-Ar(^) = 

i 

It turns out that (see Section|4]for a derivation), the limiting behavior of ttin is given by the limiting behavior 
of the unique solution of the following system: 


N 






N 




„(fc) 


(2.13) 


for k B Nr- 
We denote 


9i = 






(2.14) 


Thus, to show the self-consistent system (12.131) gives the correct limiting behavior, our goal is to estimate 
the family of quantities defined below. 

Notation 2.3. Define the 2 -dependent quantities 


Ad := max0j \Gii - , 

i 


(2.15) 


and 


Aq := max 


-y/ BiQjGi 


A := maxIAo, Aq} . 


(2.16) 

(2.17) 


We also set the control parameter 


$ := — + 


1 

■s/Nrj 


(2.18) 


The existence and uniqueness of a solution to (12.131) will be proved in Section [5l Throughtout this 
paper, we denote {u^^'^)keNr m to be the solution to (I2.13|) without specification. Our main theorem 
concerns the asymptotic behavior of eigenvalues of H in the bulk of the spectrum, that is, where the 
asymptotic distribution of the spectrum has a positive density. To make this precise, we introduce the 
following definition. 


Definition 2.4 (Bulk interval). Assume that H is a generalized Chung-Lu ensemble. Let {u^'^'^)keNr 
m to he the solution to ()2.13|) . We call a bounded interval / cz R a bulk interval if there is a positive cj 
which does not depend on N such that Im 771 ( 2 ) > c/ on {E + ip : E e 1,0 < r] 1}. 


Such an interval might not exists for an arbitrary generalized Chung-Lu ensemble, but it does exist in 
many cases of interest, for example, when (7|^^)fc6Nr,ieNN has a limit distribution, see Section|31 

Before stating our main theorem concerning the matrix model H, we need to introduce some notation. 
We follow the conventions of m- 
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Definition 2.5 (Stochastic domination). Let 

X = ■. N and Y = (y^^\u) : N e N, u e (2.19) 

be two families of random variables with Y nonnegative. We say that X is stochastieally dominated by Y, 
uniformly in u, if for all positive eonstants e and D, we have 

sup P[|X(m)| > N’^Yiu)] ^ N-° (2.20) 

ueum 


for large enough N. Moreover, we denote this by X <Y. 

Now we state our main theorem below. We point out that Theorem 1.6 of [4] contains the non-sparse 
case of our theorem, that is when the sparse parameter q = N. In fact, their theorem translated to this 
special case is much more general, since it not only deals with the bulk, but the edges as well. 

Theorem 2.6 (The local law). Let I be a bulk interval (as in Defintion \2.4^ . Fix positive d and define 
an N-dependent spectral domain 


Vi := {E + irj: Ee /, ^ rj ^ 10} . 

Then, on the domain we have A < $. 


( 2 . 21 ) 


As an application of the local law, we have bulk universality in terms of the n-point correlation functions 
of eigenvalues. The non-sparse q = N was done in Theorem 1.15 of [1]. 

Let be the n-point correlation functions of the eigenvalues of H and be the density on I. We 

(n) 

denote Pqq^ to be the n-point correlation function of eigenvalues of a GOE and psc to be the density of the 
semicircle law: _ 


Psc{E) 


V[4-ig^] + 
2 ?! 


( 2 . 22 ) 


Theorem 2.7 (Bulk universality). Let O e C'”(R”) be a test function. Let I be a bulk interval. Fix a 
parameter b = for arbitrarily small c. We have. 


lim 

TV—>cx) 


1 


ipsciE)y 


Pgoe ( + 


p{E) 
0.1 




Psc{E) ’ ■ ■ 


,...,E" + 


Np{E) 


Np{E)) 


Psc{E) J 


, , dE' 

doi .. .da„-^ = 0 


(2.23) 


for any E" G (—2, 2). 


3. Examples of interest 

In this section, we give specific examples of models that motivate the ensemble we defined in Section [21 In 
particular, we consider the cases where (7i^^)feeNr,i€NAr converges to a limit (in a sense to be made precise) 
and when the random matrix ensemble is given by the adjacency matrix of a random graph. We show that 
these models satisfy the assumptions of our theorems; note that Theorems 12.61 and 12.71 require that there 
exists a bulk interval, which is any interval where the spectral density has a lower bound independent of N. 
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(3.1) 


3.1. When ( 7 !^^) has a limit. For k e Nr, we define 


Jk 


N 


{x) :=yiri-i 


.)7i 


(fc) 


If fk in some sense, we expect the following equations from (12.131) : 


Uk = - 

mo = - 



fk{x)dx 

^ + 111=1 fk{x)uk 


2 + TIk=l fk{x)uk 


(3.2) 


where ^ 

inf fk{x)^l and/" fkix)‘^dx ^ M. (3.3) 

We note that this system of equations for fe = 1 has been derived before in [33] using the replica method. For 
(13.21) . which is the limiting, integral version of (I2.13|) . we can prove uniqueness, existence, and holomorphicity 
of the solution. 


Theorem 3.1. For any z e C"*", there exists a unique (ui,..., Ur, mg) satisfying equations (13.21) . In addition, 
the function (ui,..., Ur, mo) is holomorphic in z. 

Proof. The proof is a trivial modification of that of Theorem lA.il ■ 

Theorem 3.2. If f^ f in T^([0,1]), then any bounded closed interval I ez M. on which Immo > 0 is a 
bulk interval. 


Proof. This theorem is a consequence of the stability of the system 


(13.21) . See Theorem 15.31 ■ 


3.2. Power law with (3 > 3. For simplicity, we assume r = 1, that is, 
that 



the matrix Sy has rank one. Assume 


(3.4) 


where 0 < /x < 1/2, then the sequence of functions 


N 


f^{x) = Y,h 




(3.5) 


converges to f{x) = a;“^ in L^([0,1]). From the previous subsection, we know that bulk intervals exist. 
Therefore, Theorem 12.61 and Theorem 12.71 apply to this case. 

In particular, this model gives a random graph that has a power law as the degree distribution with 
exponent /3 = 1 + ^ > 3. To be precise, we consider a random graph with N vertices and the probability 

of the i-th vertex connects to the j-th vertex is qsij = (^) ^ ^■ Assume that different edges 

are independent. Thus, the adjacency matrix of this graph is a random matrix with variance qSij. Now 
the expected degree of the z-th edge is Jfij proportional to (^) ^. Thus the number of edges that has 

degree [Nx, N{x + da;)] is asymptotically Nx~^~'i^, up to a normalizing constant, which is exactly a power 
law distribution with exponent d = 1 + — > 3. 










4. Proof of the local law 


In this section, we prove Theorem 12.61 

4.1. Some notation and resolvent identities. We introduce some notations which will be useful later on. 
Definition 4.1 (Minors). For Tc {1,...,A^}, we define by 

:= hijliffljfj:. (4.1) 

Also, the corresponding Green’s functions are defined as 

(4.2) 

In a similar way, we use the notation 

(T) 

S - S ■ (4.3) 

i i€A(\T 

Moreover, we abbreviate ({i}) by (i) and (T u {i}) by (Ti). 

Definition 4.2 (Partial expectation). Let X = X{H) be a random variable. We define Qi by 


QfiX) := X - E [X I . 


(4.4) 


The following lemma is frequently used in the proof of the local law. See Lemma 4.5 in m for reference. 


Lemma 4.3 (Resolvent identities). For any Hermitian matrix H and T c {1,...,A^}, the following 
identities hold. Ifi,j,k ^ T and i,j i=- k, then 


r<(T)r"(T) 

^(T) ^ ^(Tfe) '^kj 

'^ij '^ij + (T) 

^kk 


and 


Gl 


(T) 


1 


G 


(T)^(T) 
ik ^ki 


^(T)^(Tfc)^(T) ■ 

'^li '^li '-'kk 


(4.5) 


If i,j ^ T satisfy i ¥= j, then 


G\ 


(T) 


(Ti) 


(Tj) 


_ ^(T) V h V cW) h 


kj 


(4.6) 


Ifi^T, then 


(Ti) 

_ 4. .. V* 4. /-i(T4)4, 

(T) “ ^ii ^ Zj ^ki ■ 

j,k 


(4.7) 
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4.2. Derivation of the limit equation. 

2 , that is, 


For the sake of simplicity, we only consider the case where S has rank 




(4.8) 


Higher ranked cases can be handled similarly. We can also recover the rank-one case by setting = /3i, for 
all i. Using the Schur complement formula (see Equation (14.7p of Lemma r4.3|) . we have 


1 



—z — 




SikGkk + Ri , 


(4.9) 


where the error term Ri is defined as 


R^ ■■= hu + Yi -Q^Y ■ 


(i) 


G 


(4.10) 


Define 




After rearranging (14.91) . we get the following equations 

1 


Ctk 


u =—y — 

N^z + akU + PkV-R^ 


1 


Pk 


V = —^y — 

Nj^z + akU + PkV-R, 


= --Y 


niN 


1 


N ^ Z + Q-kU + PkV ~ Ri 


Formally, we neglect the error terms and replace Gu with gi, then Equation (14.91) becomes 

~ ^ ^ , ^ik9k • 


(4.11) 


(4.12) 


(4.13) 


(14.121) becomes 


1 ^ 

2 


ak 


N ^^z + aku{z) + Pkv{z) 


N 


= “iv S 


Pk 


+ aku{z) + Pkv{z) ' 


1 ^ 
= -lu E 


(4.14) 


N ^^z + aku{z) + Pkv{z) 

In Section [SI we will prove that this system has a unique solution and the solution is stable under small 
perturbation. 
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4.3. Large deviation estimates. To control the error terms in the self-consistent system, we need several 
large deviation estimates, which we state here. We omit the proof here, which is a slight modification of 
Lemma A.l in m- 

Lemma 4.4. Assume that the family of random variables Equation (12.61) . 

(i) Fix an index i G Nat. Let {Ai,..., An) be a family of random variables that is independent of 
(hii ,..., hiN) and satisfies Aj < 1 for j G Nat, then we have 


'^{hij{t) Sij)Aj 
3 



(4.15) 


(ii) Fix distinct indices i,j G Nat. Let be a family of random variables that is independent of 

{hii ,..., hiN) ciiT'd {hij ,..., hNj) and that satisfies Bf^i < 1 for k,l e Nn, then we have 


X! hik{t)Bkihij{t) < \j~^ + \Bki\^'\ ■ (4.16) 


4.4. Proof of the local law. Our strategy is to first assume that A is bounded by a large control parameter 
N~‘^ on some event that has high probability. With this a priori bound on A, we prove that A is actually 
bounded by a much smaller control parameter $ (see Notation 12.31) with very high probability. Thus, the 
probability that A lies in between N~‘^ and 4> is very small. However, when r] is large (order one), we can 
easily show that A is bounded by <i>. Finally, we use a continuity argument to push this estimate all the way 
down to ?7 > for arbitrarily small positive 5. 

The following lemma essentially says that with very high probability, A is outside the interval [$, 

So [4), 7V““] is sometimes called the “forbidden” region. 

Lemma 4.5. Assume that L is a bulk interval (see Definition \2.4\ ) ■ Let (f be the indicator function of some 
(possibly z-dependent) event. For any spectral domain T> c Bi, if tfA < for some positive c, then 

(j)A < 4). 

Proof. On the domain 2?/, it is easy to see that u and v are bounded and Imn and Imu are bounded 
below. By the definition of A and V^, we see that 


where 9i is defined as in (12.91) : 

First, we bound Aq. We need the identity 


4>Gii 


ai + Pi. 


(b) 




(4.17) 


(4.18) 


(4.19) 
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which is easily obtained by iterating Equation (14.61) of Lemma 14.31 . We get an a priori bound on using 
Equation (14.51) of Lemma 14.31 


< 4‘Gjj + 


) _ ^ 


Similarly, one can prove that 


By (ii) of Lemma 14.41 (j4.191) yields 




Uit/j 


19i9j 


+ 


I9i9i „ „ "tiij) ^ 






k,l 


G 


kl 


By the Cauchy-Schwartz inequality and Ward’s identity, we have 


>Y,^k9i 


k,l 


G 


(b) 

kl 


< 




92lmG^f^ N 


< — 
1 


(4.20) 


(4.21) 


(4.22) 


(4.23) 


where we used the bounds (l>G^kk ^ ^ Sfe ^ Therefore, Equation (14.221) is stochastically 
dominated by 


J_ f /Ml + /Ml I = M_$. 


(4.24) 


So Ao < $. 

Next, we bound the error term Ri. First, ha < 1/s/q < 4>. Second, using Equation (14.171) and the bound 
Ao < 4>, we see 


, v"! GikGki . v-i 
92_jSxk „ = 92_j^ik 


(i) 


G., 


GikGki 

G,. 


+ (bsiiGii < — y — + - < 9i^ . 
“ TV V 6»i q 
k 


k “ k 

Third, we use Lemma [4.41 and the Ward identity (as we did in (14.231) ') to estimate 

(i) i^) (i) 

<t>Q^ 2 h,kG^^}hu = cl>Y,{hl - s.k)Gl:l + <). ^ h,kG^^^hu 

k,l k k^l 


< ^ + 


9, 


Vg V^Vn 




(4.25) 


(4.26) 


Therefore, (j)Ri < 9i^, which yields 


Gii + 


1 


Z OliU -\- PiV — Ri 


< $ . 


(4.27) 


By the stability of the self consistent equation. Theorem 15.31 we have cj)9i \Gii — gi\ < 4) and therefore, 
A < $. ■ 
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For r] large, we establish an a priori bound on A, which we use as an input for Lemma 14.51 


Lemma 4.6. For any bounded interval 1 1 


we have A < $ on the line segment {E + ir] : E e I,r] = 10}. 


Proof. We proceed as in Lemma HTSl First, we want to show that Gn < In Lemma HTSl this estimate was 
obtained using the a priori assumption A < but now we only have the trivial bound \Gii \ ^ ^ = 0(1). 
First, note that 

(i) 

k^l 


< 


EL 


(4.28) 


Next, we proceed with Equation (14.7|) . to see 


\Gu\ = 


1 

^ ^ + (^) ■ 

Next, we get a lower bound for ImG^.^- Using Lemma 03] and Equation (I4.17L we see 




V + hli ImG^if + 0< ^ ry + 0< 




O<(0f) 


(4.29) 


(4.30) 


Thus, ^ < ImG^E Note that the same holds for Gkk. Now, we return to Equation (I4.29L to get |Gii| < 

Once we have this bound, we can move on to estimate G^ < as we did in ()4.22p . The estimate 

for Ri proceeds in exactly the same way as in Lemma [4.51 we omit the details: Ri < 9i^. Thus, 


Gii + 


1 


Rr 


$ 




z -\- (y.iU -\-13iV^ — Ri 

By Theorem I5.5[ this implies 6i \Gii — gi\ < which in turn implies < $. Therefore, A < $. 


(4.31) 


Now, we complete the proof of Theorem 12.61 with a continuity argument. 

Proof of theorem 12.61 We choose a lattice A c Dg such that |A| ^ and the A“®-neighborhood of 
A covers the whole Dg. Take <() = [A ^ for some positive and small enough constant c. Lemma [4.51 

states that for any positive, large D, e small (such that < N~G, and N large enough. 


¥[3w G A : ^ A sS N-^] ^ N-° . 

(4.32) 

Since A is Lipschitz on Dg with constant at most we have 


F[3z G A : 2$A® ^ A ^ A"^] ^ A"^ . 

(4.33) 

Using Lemma [4.61 and the continuity of A, we have 


Ppz G A : A 5= 5>A^] sg N-^ . 

(4.34) 

Thus, A < 4). 

■ 
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5. Solution to the self-consistent equation 


In this section, we analyze the self-consistent system (EH). We prove the solution exists and is unique, we 
prove the equation is stable, and we prove some asymptotic behavior of the solution. 

For the sake of simplicity, we only consider the two dimensional case, as higher dimensional cases can be 
handled similarly. Thus, we have the equations 


u{z) 


v{z) 


m{z) 


1 

N 




_Ofe_ 

2 -t aku{z) + I3kv{z) 


1 y /3fc 

N ^^z + aku{z) + Pkv{z) ' 


1 

N 




1 

2 + afcu(z) -I- l3kv{z) 


(5.1) 


This is a two-dimensional version of equations (12.131) . Next, we have an existence and uniqueness theorem 
for the solution of equation (EH). 

Theorem 5.1. For any z e C"*" there is unique (u,v) G C"*" x C"*" satisfying (EH- Moreover, (u,v) is 
holomorphic in z. 

Remark 5.2. Thus, m is uniquely determined by (it,u), by the third equation of (15.11) . 

Proof. Fix z g C+. Define F : C+ x C+ C+ x C+ by 




/1 wi ckfc ^ y* 

yN ^ z + ako- + PkT' N ^ z + aku + PkT 


(5.2) 


As a rational function, it is continuous. By Brouwer’s fixed point theorem, there is a fixed point (u, v) G C"*" x C"*" 
such that F{u,v) = {u,v). 

To see uniqueness, we first note that by taking the imaginary part of EH that 


1 afe(afeImM-I-/3felmu) 

Imu ^ ^— 

N , U 4- rvi.7/. -I- 


finrl 


1 Bh(ah Im u + Bh Im v) 


(5.3) 


For T defined by 


T := — y 

N ^ 


k \z + otku + I3kv\ 


^k 


Pi _ 


(5.4) 


the Perron-Frobenius theorem implies the spectral radius r(T) is less than or equal to 1. 

Now, assume for the sake of contradiction that there is another solution {u',v'), then again we have an 
analogue of Equation (15.41) and r{T') ^ 1, where T' is defined likewise: 


So, we have 


T'-.= -y 

N ^ 


^ \z + akU + Pkvf 


al akPk 
O-kPk Pi 


= ly 

N ^ 


ak{ak{u - u') + j3k{v - v')) 


— - 

N ^ {z + Uku -I- I3kv){z -I- aku' + (Bkv') 


(5.5) 


(5.6) 
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and 


Denote 


/ 1 v* Pkiakiu - u') + I3k{v - v')) 


= -y 

N /-J 


N ^ {z + akU + Pkv){z + aku' + Pkv') ' 


1 


f := — y 

N ^ 


1 


N ^\z + akU + f3kv\ \z + Uku' + j5kv'\ 


ai akPk 
Oikhk Pi 


and ^ = (|m — ti'l, I?; — f'l). We take the absolute values of equations (15.61) and (15.71) . thus 

^ = n- 

Then we take the Euclidean norm to see 


lien < 


Ti 


However, applying Cauchy-Shwarz inequality to dSH) yields 


n 


< 2(reii + ||n||)<iieii 


(5.7) 

(5.8) 

(5.9) 

(5.10) 

(5.11) 


which is a contradiction. 


We also have a stability theorem in the bulk. This theorem says that the solution to equation (15.11) is 
stable under small perturbation of the equation. 

Theorem 5.3. Let {u,v) be the solution to (EB. Assume / G R is a bounded internal and c > 0 is a constant 
such that Im m ^ c on 

Vi = {z :RezB I ,lmze (0,10)}. (5.12) 

Then there are e,C e R"*", depending on c but not on N, such that the following holds. For any z G Vj, 
assume that {u',v') satisfies 

1 

- =ri 

) 


u 


N ^ z + aku' + Pkv' 




Pk 


(5.13) 


7 =^2 


N ^ z + aku' + PkV‘ 
and max{|ri|, |r 2 |, |u' — m| , |i;' — u|} ^ e, then 

max{|u' — it| , — ti|} ^ Cmax{|ri|, |r 2 |} . 

Proof. Denote 


f := — y ■ 

N ^ 


N ^ {z + akU + Pkv){z + Uku' + Pkv') 
Then, subtracting (15.2211 from (15.IL 

{I -f){u- u\v - f) = (r-i,r 2 ). 


al akPk 

OikPk Pi 


(5.14) 

(5.15) 

(5.16) 


The assumption 1mm > c implies that rt, u', v, v' are bounded and Imu, Imu', Imz;, Imw' are bounded below. 
Denote 




1 




1 


^k ^kPk 


1 V'l 
























Then 


Ti-T 


^ C 2 max{|w' — m| , \v' — Thus we have 


max{|M' — u\, \v' — u|} < \\{I — 2\) ^|| (max{|ri|, |r 2 |} — (max{|u' — u\ , \v' — 'y|})^) ■ 


(5.18) 


It remains to show that the eigenvalues of Ti are bounded away from 1, so that {I — T) ^ is uniformly 
bounded. This, however, is guaranteed by Lemma 15.41 ■ 


Lemma 5.4. K c C"*" is a compact set. Assume that A = 
the spectral radius of 

lAI := 


ail ai2 
021 0,22 


, where Oij e K,i,j e {1,2} and that 


(5.19) 


|aii| |ai 2 | 

|a2i| |a22| 

is ^ 1. Then there is a S > 0 depending only on K such that the eigenvalues of A are bounded away from 1 
by 5. 


Proof. Assume that this is not true, then we can take a sequence A„ satisfying the same conditions as A, 
with A„ and f,n A 0 such that 

An^n = with A„ ^ 1. (5.20) 

We take a subsequence to get A and ^ such that 


= (5.21) 

By Perron-Frobenius theorem, we know that \A\ |^| = |^|, which means On^i is parallel to ^i. This contradicts 
with Oil G K. I 


We also need a stability theorem when rj is large. The difference is that we do not assume an a priori 
bound for max{|ri|, |r 2 |} here. 

Theorem 5.5. Let (m,u) be the solution to (EU. Assume / G K. is a bounded internal. Then there are e,C e 
K"*", depending on c, / but not on N, such that the following holds. For any z G {() : ReC G /,Im^ = 10}, 
assume that {u',v') satisfies 



_Ofc_ 

2 ; + aku' + Pkv' 

Pk 

z + aku' + j3kv' 


ri 


r2 


and maxljw' — u|, \v' — u|} ^ e, then 

max{ \u — u I , |i;' - u|} ^ Cmaxllril, |r2 |} . 
Proof. As in the last theorem, we get 

{I - f){u - u',v - v') = (ri,r 2 ). 

Recall the definition (15.51) of T'. We take the imaginary part of (15.221) to get 


Im u' 



r)Oik 

+ Imri 


\z+aku'+Pkv'[‘‘ 

Im v' 


Us 

vPk 

-1- Imri 


\z + akii' 


(5.22) 


(5.23) 

(5.24) 


(5.25) 


When ri and r 2 are small enough, the right hand side is positive, thus r(T') ^ 1 by Perron-Frobenius 
theorem. We already know r{T) ^ 1, thus r(|T|) ^ 1 (see (I5.19P for the definition of |T|). By Lemma [5.41 
the eigenvalues of T are bounded away from 1, thus {I — T)~^ is uniformly bounded. I 


16 






















6 . Bulk universality 


In this section, we prove the universality of n-point correlation functions. The key ingredient is a universality 
theorem for deformed GOE matrices [3T], which we state below. 

A matrix V is said to be (I, G)-regular at E if for some parameter I such that 

— sS / ^ and iV'^U ^ ^ (6.1) 

there exist positive constants cy and Gy such that 


Cy ^ ImTOy(i?' + ir]) ^ Gy (6.2) 

uniformly for E' e {E — G,E + G) and I ^ rj ^ 10. In our case, V is H diagonalized and I = so 

T = for arbitrarily small e' > e > 0. Here is their theorem. 

Theorem 6.1 (B. Landon & H.-T. Yau [31]). Let 

H=V + Vfw, (6.3) 


where V is a (deterministic or random) diagonal matrix and W is a standard GOE matrix. Suppose that V 
is {I, G)-regular at E and that N~^G‘^ > T > N^l for some positive e. Let O e G®(IR.”) be a test function. 
Fix a parameter b = for any positive c satisfying c < e/2. We have, 


pE-\-b p 

/ / 0 (ai,...,a„) 

J E-h 


pE+h 

lim 




E' + 


ai 


^ AO I 7-.// , Q!i 


ipsc{E)y 




AN) 


m' 


E" + 


Np%\{E) ' 

(Xin 


-,...,E' + 




pic\E)^ 


dai ...dun^^ = 0 . 


(6.4) 


Our strategy to prove bulk universality is as follows. Let Ht = {hij{t)) be a symmetric N x N matrix. 
The dynamics of the matrix entries are given by the stochastic differential equations 


dhij{t) = 


dBjj (t) _ hij{t)dt 
VN ‘2NSrj 


(6.5) 


where B is symmetric with {Bij)i^i^j^]\f and {Bu/V2)(L^ a family of independent Brownian motions. The 
initial data Hq is our target matrix H. We shall prove a comparison theorem, Theorem 16.61 which states 
that, for T = with e « 1, the spectral statistics of Ht agree with those of Hq. Next, we show that 

Ht can be written as the sum of two random matrices, 

Ht = HE-dG, (6.6) 

where Ht has a local law iTheorem 16.71) and dG is a small GOE matrix. Next, we apply Theorem 16.11 to 
get bulk universality for Ht. Thus Theorem 16.61 savs that universality also holds for Hq. 

We start with the comparison between Hq and Ht. The following lemma is a slight modification of 
Lemma A.l in [23) . The main difference here is that we do not assume Sij ~ 
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Lemma 6.2. Assume Sij are as in (lO) . Suppose that F is a smooth function of the matrix elements {hij)i^j 
satisfying 

sup E ((iVs-.i < M, (6.7) 

where (9H)ij = Oijhij with O^i = 1 unless {fc, Z} = {i,j} and 0 ^ 6ij ^ 1 , then 

EF{Ht) - EF{Ho) = 0(tN)M. (6.8) 

Remark 6.3. In our case, t ^ thus we only need M < for some c ^ £. 

In order to apply this lemma to our case, we need some uniform bound on the Green’s functions 
Gt = {H — z)~^ for t G [0,T]. Clearly, for every t > 0, we have E{hij{t)^) = stj. Moreover, hij{t) is a 
small perturbation of the initial hij ( 0 ) and has the moment bound 

Thus, we expect the same local law to hold for Fix- 

Theorem 6.4 (The local law for Ht). Assume that I is a bulk interval. Fix positive S and define an 
N-dependent spectral domain := |i5 + iTy : i? G ^ rj ^ 10}. Let {hij{t))^j^^ be defined as above 

in (16.5|) and fix t G [0, Therefore, on the domain Fg, we have A < <i>. 

The following comparison lemma is a consequence of Lemma 16.21 and Theorem 16.41 It is a modihcation 
of Lemma 5.2 in |30) . 

Lemma 6.5. Assume that I is a bulk interval. Let Ht be defined as in dSH). Let positive 8 be arbitrary and 
choose an rj such that N~^~^ ^ rj ^ N~^. For any sequence of complex parameters Zj = Ej±iri,j = 1,..., n, 
with \Ej\ G I with an arbitrary choice of the signs, we have the following. Let Gt{z) = {Ht — z:)~^ be the 
resolvent and let F{xi,X 2 ,. •., Xn) be a test function such that for any multi-index a = (oi,..., a„) with 
1 ^ |a| ^ 3 and for any positive, sufficiently small uj, we have 


max||d“F(a;i,...,a;„)| :max|xj| < ^ (6.10) 

and 

max I \d°'F{xi, ..., x„)| : max \xj \ ^ ^ N^° (6-11) 

for some constant Cq. Therefore, for any 8 with N~^~^ ^ 77 ^ N~^ and for any choices of the signs in the 

imaginary part of Zj, we have 


E[F 


(4TrG,(.,) 


1 

N 


TrGt{zn) 


- E[F(Gt 


Go)]] 




GtN^+‘^^ 


where c and G are constants depending on Cq. 

Proof. We consider only the n = 1 case for simplicity. We want to show that 


N 


E[F -TrGt(zi) -E[F -TrGo(zi) ]] 


N 


< 


CtN^+^^ 


( 6 . 12 ) 


(6.13) 
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We need to compute the derivative of i TrG with respect to the (i, j)-th entry. Note that 


5a\^TrG 




(6.14) 


where J is the matrix whose entries all vanish except the (a, b) and ( 6 , a)-th entries. The right hand side is 
a sum of products of off-diagonal entries of G and the number of terms is order N. Thus, we need a bound 
on the off-diagonal entries of G, down to the scale rj > N~^~^. In order to get this bound, we first derive a 
delocalization bound for the eigenvectors in the bulk. For Xk G I, we set E = Afc 


Im Gii = ^ 


I 7 |2 

l<l V 

|Aa - zf 


> 


(6.15) 


By Theorem 16.41 ImG^^ is bounded uniformly in on an event with probability 1 — N~^. Hence on this 
event, we have as long as Afc e /. This still holds if we replace / with a slightly bigger interval 

I (independent of N), because of the analyticity of m near I. Now we estimate, for Eel, that 

1 


7 h 


7 h 




-- V - 


-+ V - 



^ CN-^+^ 2 

i Xisi XifI Xiel 

We define a dyadic decomposition of the eigenvalues lying in I: let 

Uo := {j : A, e /, |A, - E\ ^ N-^+^} , 

C/„ := {j : Xj G < \Xj - E\ ^ 2^N-^+^} , 


|Ai - z| 


+ C. 


(6.16) 

(6.17) 

(6.18) 


for 1 ^ n sg log 2 N, and 

t/iog, 7 V +1 := {j : A, ei,N^ < \X, - E\}. (6.19) 

With probability 1 — N~^, we have |G„| ^ , for 0 ^ n ^ log 2 N, which can be seen below. For any 

J e i such that | J| > , we set E' to be the middle point of J and 77 = then we have 


1 


G > Im m{E + irj) = — ^ 


N^iXi-E^ + r, 

The dyadic decomposition enables us to estimate 

\gme + iv)\^cn-^+^Y. 2 

n ieUn 


> — 

'2 IV 


XieJ 




|J|iV 


#{J : A, e J }. 


\Xi- E - iri\ 


^ GN' 


35 


( 6 . 20 ) 


for rj ^ N ^ ^ and therefore, ^ |Tr(9^{,G| ^ on an event with probability 1 — ^. Outside this 

event, we have a deterministic bound |Trd^f,G| ^ G7V2(i+^). apply Lemma [6.21 to complete 

the proof. ■ 

Lemma 16.51 readily implies the following comparison theorem for correlation functions: 
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Theorem 6.6. Assume that I is a bulk interval. Let Ht be defined as in (16.51) and T = N with e 
arbitrarily small. Let and be the n-point eorrelation functions of the eigenvalues of H andHx, then 
for any test function O e we have 


lim [ 0(ai,...,a„) p^fP^ (e + 


N'' 


,E + 


N 






)] 


= 0 . 


( 6 . 21 ) 


Proof. This theorem is a consequence of Theorem 2.1 in [23] and Lemma 16751 . 


It remains to show that Hr has bulk universality. We note that since Sij > c/fV, Ht has a small Gaussian 
component, that is, Ht = H + cy/TG. We state without proof the local law for H here: 

Theorem 6.7 (The local law for II). Assume that / cz R is a bulk interval. Fix positive S and define 
an N-dependent spectral domain Hg := {E ir/ : E e I, ^ p ^ 10}. Let {hif) be defined as above. 

Therefore, on the domain Hg, we have A < $. 

Proof of Theorem 12.71 By Theorem 16.11 and Theorem 16.71 the conclusion of Theorem 12.71 holds for Ht 
in place of H. Then, by Theorem 16.61 the same conclusion holds for H. ■ 


7. Adjacency matrix without random signs 


In this section, we consider a different matrix ensemble derived from the adjacency matrix of a random graph 
without random signs (recall Equation (12.2^ 1 and outline the proof of a local law and bulk universality. We 
consider A = (oy) such that 


P [fly = 1] = ^ and P[ay = 0] = 1 - ^ , 


(7.1) 


where Sij has the form <\Tn . Note that Eoy A 0, so we centralize A by subtracting its expectation EA = 
normalize by We denote the centered, normalized entries by /ly. Thus, /ly has law 


hij — 


■s/q ^ 




N ) 


qsjj 

N 


and P[/iy = — 


1 gsy 

Vq N 


= 1 - 


qs^ 

N 


(7.2) 


As before, using the Schur complement formula (I4.7|l . we find the equation SikGkk + 

where the error term Ri is defined as 


(d 


R^-.= hu+Y, E • 


(7.3) 


Thus, we get the same self-consistent equation for u, neglecting the error term. We need to bound the 
additional error term s^k^kk- We have, on the event (j)K < N~^ (see Lemma H3]) . that 


R^ < ^ T 



N \ N 


< 


(7.4) 


If 9 ^ VN, then Ri < -^ < Oi^. Once we have this bound, together with large deviation bounds as in 


Lemma l4.41 we readily get the local law, repeating the proof of Theorem l2.6l In summary, we have Theorem 

rm 
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Theorem 7.1. Assume that q = with k ^ 1/2. Let hij be defined as in (EH). Assume that / c R is a bulk 
interval. Fix positive 6 and define an N-dependent spectral domain := \^E -\-irj : E e ^ rj ^ lO}, 

then on this domain A < $. 


As in the original model, the universality in the bulk holds for this model. We state the universality 
theorem without proof, since the proof for Theorem 12.71 can be reproduced almost verbatim. 

Theorem 7.2. Assume q = N'^ with k e (0,|]. Assume that / c M is a bulk interval. Let be the 

n-point correlation functions of the eigenvalues of F[ and be the density on I. Let O G be a 

test function. Fix a parameter b = for arbitrarily small c. Therefore, we have 


r ^ / /-\i 

L-b L ■ ■ ■ ’ V ^ ■ ■ ■ ’ ^ 


^E+b 

lim 

JE-b 


A 

We)) 


iPsciE)) 


^ ^(”) I pii I 

Pgoe E + 


pic\Ey 


,...,E" + 


pic\E)^ 


, , dE' 

dai . ..dan-^ = 0 


(7.5) 


for any E" G (—2, 2). 

Remark 7.3. With some more effort, we can remove the restriction k ^ 1/2 in Theorem 17. II and Theorem 
m The idea is to solve a general self-consistent equation 


9i = 


1 




N ^ 
k=l 


(7.6) 


for i G Nat, where Sik = Sifc(l — ^|p). We refrain from proving this in any detail here, as the most interesting 
case in applications is when k is small, that is, when the graph is very sparse. 


A. General self-consistent systems 


In this section, we forget about N and consider the general self-consistent equation 


9z{x) 


1 

^ + {S9z){x) 


(A.l) 


for 2 G C+. Eouation (lA.ll) has earlier been studied in 0, m and [53] . Recently, it was extensively studied 
in 0, and named quadratic vector equation (qve). We shall prove a theorem about the uniqueness and 
existence for the solution. Our theorem is contained in much more general results in 0 and [29] . For the 
readers’ convenience, we present here a different but short proof in a much simpler setting. 

Denote m{z) = g(x)dx. We let x G [0,1], g^ G L^([0,1], C+) and Sg{x) := fg s(x, y)g(y)dy with 
s G T^([0,1]^), s ^ 0 on [0,1]^. 

Theorem A.l (Existence and uniqueness). For any fixed z g C+, there uniquely exists a Lf function 
gz '. [0,1] —> C"*" that satisfies equation (jA.lll . Moreover, g as a function from C"*" —> L^[0, 1] is holomorphic. 

Proof. Existence. Our strategy is, we first find a sequence of functions gk satisfying the equation at dyadic 
points X = ^,0 I 2^. Second, we show that the sequence is bounded in L^[0, 1]. Third, we use the 
fact that S' is a compact operator to find a convergent subsequence. Fourth, the limit turns out to be the 
solution. 
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Define a map : (C+)2*’ —> as follows: For a vector ^ ^ define a staircase 

function 

2 *’ 

^c(2;) = 2 1[L:^,4](2;)C' ■ (A.2) 

1=1 "" "" 

Define a projection operator Pk on D^[0,1] such that for v G D^[0,1], 



v(t)dt, X G 


i-i r 

2 ^ ’ ^ 


(A.3) 


The projection Pk is actually an operator from 




So we define 


The function (j)z is continuous on cjiz '■ (C+)^*’ —> and |(/) 2 |go is bounded by i. By Brouwer’s fixed 

point theorem, there exists a C* G C+ such that (pziC*) = C* and that |C*loo ^ 

We set gk ■= By (IA.4I) . we have gk = Pk ^ z+Sg ^ )' function gk is bounded by i, hence is 

bounded in T^[0,1]. We can find a subsequence, still denoted by gk, such that Sgk g in L^. Note that 


-1 


9k 


z + Q 




-1 


-1 


z + Sgk z + Q 


< c\\9k-g\\ + 


{Pk - I) 


+ {Pk - I) 

-1 


Z + Q 


Z + Q 

0 . 


(A.5) 


Thus, gk converges to g := , which satisfies g = — ■ 

Uniqueness. Fix z = E + ig e C+. Assume that there are two different solutions g and g. Denote 


a(x) = z + (Sg)(x) (A.6) 

and d similarly. Since Im 5 (a:) > ; we have 

|a(x)| 

la(x)llmg(x) > f la{y)llmg(y)dy. (A.7) 

Jo \a{x)a[y)\ 

Next, we use a generalization of Perron-Frobinius theorem. 

Theorem A.2 (The Krein-Rutman theorem). Let X be a Banach space, K X a total cone (that is, 
K — K = X) and T : A —> A a compact linear operator that is positive (that is, T{K) cz K) with positive 
spectral radius r((T), then r((T) is an eigenvalue of T with an eigenvector u G A\{0} : Tu = r(T)u. 

In our case, A = L\0, 1], A = {/ G A, / ^ 0}, and T is the integral operator with kernel The 

operator T is compact because its kernel is continuous. The theorem tells us that there is a u{x) G A\{0} 
such that 

/ " r{T)u{x). (A.8) 

Jo Hx)a[y)\ 
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We take the inner product of (IA.7I) with the eigenvalue u to get <!u. a Im > riTYu. a Im q\ Thus, r{T) < 1. 


Since T is symmetric, this means ||T|| < 1. Similarly 
On the other hand 


T 


< 1 . 


\g{x)-g{x)\ = 


f s{x,y) \g{y) - g{y)\dy 
a{x)a{x) 


(A.9) 


Denote w(x) = ya(x)&(^ \g{x) — 


'j{x) sS [ 
Jo 


s{x,y) 


So Ikll 


0 ^^a{x)a{x)^ya{y)a{y) 
(T + T) ||r(;|| < \\w\\, a contradiction. 


wiy)dy ^ ( -(T + T)w ) (x). 


(A.IO) 


The following theorems concern the stability of the equation. Here we assume s to be bounded below, so 
that when 77 gets small, we still have a bound on gi. We omit the proof here, because the proof is essentially 
a reinstatement of the uniqueness of solutions. Our theorem are contained in more general theorems in 
We present our theorem here for the readers’ convenience, since they are easier to read. 


Theorem A.3 (Stability in the bulk). Assume that s is bounded below on [0,1]^. Assume g is the 
solution to (lA.lll . Denote miz] = I)? Q(x)dx. Let I be a bounded interval such that Imm(z) is bounded below 
on V := {E + ir] : E e 1,0 < g ^ 10}. 

Thus, there are positive constants C and Eq such that for any z B D, if g satisfies 


g + 


1 


z+ Sg 


^e<£o,||5-5z|| <eo, 


(A.ll) 


then on V we have II 5 — gz\\ ^ Cs. 


Theorem A. 4 (Stability in S). Assume that s and s are bounded below on [0,1]^. Assume g is the 
solution to (|A.1I) and g solves (|A.1I) with S replacing S. Let I be a bounded interval such that Imm is 
bounded below on D := {E + ig : E B 1,0 < g ^ 10}. Thus, there are positive C and Eq such that for any 


z B D, if S satisfies 


S-S 


^ E ^ Eq, then on V, we have \\g — g\\ ^ Ce. 


B. Local law and bulk universality of a general ensemble 


In this section we consider Sij e [c/iV, C /for some c,C > 0, and do not put low-rank conditions on (s^). 
We have Theorem IB.II and IB.31 which are parallel to Theorem 12.61 and Theorem 12.71 Again, we point out 
that in the special case q = N, Theorem IB. II and IB. 31 are contained in Theorem 1.6 and Theorem 1.15 in [4]. 
We set a sparse parameter q = N'^ with k e (0,1). Assume either the biased case where 


P 



1 " 


qsij and P [hij = 0] = 1 — qSij 


(B.l) 
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after centralization or the unbiased one where 


hij = ± 






and P \hij = 0] = 1 — qsij 


(B.2) 


In either case, we have moment bounds 

(B.3) 

So, we have large deviation bounds, assuming \Aj\ < 1 and \Bm \ < 1: 

/ \ 

(B,4) 

and 

. (B.5) 

fc#i v9 k,i ) 

As before, we set a control parameter $ := ^ + We think of as a step function on [0,1]^, 

that is, y) = Ns\^\ for x e (^, -j^] and y e Let be the integral operator with kernel 

Thus, is uniformly bounded in the Hilbert-Shmidt norm. Consider the self-consistent equation 


g^^\x) 


1 

Z+SWg{N)^^) ’ 


(B.6) 


where g : C+ x [0,1] ^ C+. By Theorem lA.li g^^'> exists and is unique. It is easy to see that is also a 
step function on [0,1], whose value on (■^, ;^] we denote by gi. Define := g^^^dx. In the rest of 

this paper, we will omit the IV for simplicity. As before, we are going to analyze the Green’s function defined 
by G := {H — z)~^. First, we apply Schur’s complement formula (14.71) . to get G~^ = —z — SikGkk + Ri, 
where the error term Ri is defined as 


R^ ■■= (B.7) 

k k,l 

Neglecting the error term Ri, we get equation (IB.61) . 

We define bulk interval as in Definition 12.41 We have a local law as follows. 

Theorem B.l (Local law for the sparse model). Let H be a family of random matrices defined 
above. Assume that ^ [ AT > ■^] for some positive constants c and G, that g = g^^'> solves the self-consistent 
equation D, and that m{z) = m^^\z) = g^^^x)dx. Assume that I is a bulk interval. Thus, on the 

domain 

Vl-.= {E + \y.Ee . (B.8) 

we have A < $. 
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Sketch of proof. Assume that (j)K < N where (j) is some event. We shall show A < <I) := . In the 

following estimates, we omit universal constants that are independent of N. First, we get some information 
from the fact that Imm is bounded below, that is, 


\9^\ ^ 


1 


1 1 

^ 


Imm ci 


(B.9) 


So g is uniformly bounded in A”[0,1]. Moreover, 

Sj Srj luvgj 


Iragi ^ 




^ Im TO ^ C/ . 


(B.IO) 


So Imgi is uniformly bounded below. Therefore \Gii\ — 1. Moreover, < 1 by resolvent identities. 
For Aq, we have Gij < <i>. For A,;, we firstly have < $, which implies 


max 


Gii + 


1 


< < 1 ). 


2 + Xij 

Theorem I A. 31 tells us that \Gii — gi^ < $. By the L®-boundedness of s, we have 

< $ . 


(B.ll) 


1 


z + SijGjj z + {Sg)i 


Therefore, 


Arf ^ max 


Gii + 


1 


z + SijGjj 


+ max 


1 


1 


z + SijGjj z + {Sg)i 


< $ . 


(B.12) 


(B.13) 


Finally, we apply a continuity argument to conclude the proof. 


Note that the theorem does not assume any limit of Instead, it requires to have some common 
“bulk,” that is, a uniform lower bound for Imm near some interval / c R. In particular, if is close 
enough to some certain s*, whose ImTO* is bounded below on T>q, then we have uniform lower bound for 

ImTO. 


Corollary B.2. Assume the same conditions as in Theorem \B.l[ We assume s* satisfies the conditions of 
Lemma \A.3l and that ImTO*(z) is bounded below on some spectral domain Dq := {E + ii] : E e 1,7] e (0, 10]}. 
For any positive 6, define Dg := {E + ii] : E e 1,7] e 10]}. Thus, there exists some positive Eq such 

that if ~ '^*11 ^ ^ 0 ; we have A < $. 

Theorem B.3. Let be a family of random matrices defined above. Assume that Sij e ■^] for some 

positive constants c and G, that g = g^^'^ solves Equation (IB. 6 |) . and that m{z) = m^^\z) = g^^^(x)dx. 

Let L be a bulk interval. Let be the n-point correlation functions of the eigenvalues of H and be 

the density on I. Let O e be a test function. Fix a parameter b = N'^~^ for arbitrarily small c. We 

have, 


pE-\-b p V ^ / rv 

lim / / 0(ai,...,a„) —— (e'+ j ...,E' + 




1 


{Psc{E))' 


:.p'S6e E" + 


pi’^\E) 


,...,E" + 


pic\E)^ 


dai . ..dan‘^ = 0 


(B.14) 
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for any E" e (—2, 2). 


The proof is parallel to that of Theorem 12.71 
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